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(54) Method for detection of elements of interest in a digital radiographic image 



(57) A method of improved detection of elements of 
interest by means of a detection system in a digital radi- 
ographic image of an object, acquired on an acquisition 
chain, comprises a first calibration phase in which the 
performances of the detection system are established. 
For such purpose, a calibration curve with several 
dimensions is determined by means of a mathematical 
model of the acquisition chain and object, in order to 
express a false positive probability as a function of a set 
of parameters corresponding to the acquisition chain 
and object. In the course of a second phase of use. a 
false positive probability value is determined for each 
pixel of the image. One then selects among the pixels of 
the digital radiographic image acquired those whose 
probability value obtained satisfies a predetermined cri- 
terion. 
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Description 

[0001] The invention concerns the automatic detection of small-sized elements ol interest mingled with noise in a 
digital radiographic image. 

[0002] It has a particularly interesting application in the medical field, in which analysis of the fixation of calcium 
salts in the tissues of the organism is undertaken. These micro-calcifications are of small size, and it is often difficult to 
distinguish them from the noise in an image, where microcalcifications and noise are represented at the same time in 
the form of small spots. 

[0003] The effective treatment of a breast cancer necessitates a medical follow-up from the first signs of cancer. 
These signs are present in the form of calcium salts fixed on the tissues of the mammary gland. They are generally 
detected on an X-ray of the mammary gland, in the course of which they appear on an image in the form of minuscule 
bright spots (light spots). 

[0004] Their small dimensions render their detection difficult and they are sometimes missed. 
[0005] Entirely digitized mammography acquisition chains are known, intended for the detection of microcalcifica- 
tions. Said detection necessitates the development of image processing algorithms in order to differentiate the light 
spots due to said microcalcifications from those due to noise. 

[0006] In general, in the prior art, the detection of microcalcifications is carried out from a film that is digitized. Fil- 
tering can then be carried out, followed by application of a threshold on the image. 

[0007] One can cite the method of estimation of noise level from the image itself, developed by Nico Karssemeijer, 
"Adaptative Noise Equalization and Recognition of Microcalcrfication Clusters in Mammograms," in "State of the Art in 
Digital Mammography Image Analysis," published by K.W. Bowyer, S. Astley, World Scientific, 1994. But this method 
leads to an ambiguity of distinction between the elements representing noise and the elements representing microcal- 
cifications, for it is based on scanned film images and does not take the acquisition parameters into account. 
[0008] Another method was presented by Michael Brady, Ralph Highman and Brasil Shepstone, "A representation 
for mammographic image processing," Medical Image Analysis. Vol. 1, No. 1, 1996, pp. 1-18, in which a process of 
image generation on X-ray film is modeled. But this process does not include a noise model. 

[0009] In all the systems of the prior art, the selection of a threshold is often difficult, since the noise is dependent 
on the image content. The size of the microcalcifications not being limited by a minimum value, their detection results 
from a compromise between the detection of many noises as microcalcifications and a neglect of small-sized micro-cal- 
cifications. 

[001 0] The invention is directed to a solution to the problem of determination of the threshold in order to extract the 
microcalcifications from the noise. 

[0011] In an embodiment of the invention an estimate is made of the probability for a pixel or an elementary zone 
to be a false negative, that is, the probability for a pixel representing a theoretical microcalcif ication not to be detected, 
in order to afford the user the choice of an acceptable number of false negatives. 

[0012] In an embodiment of the invention a method of improved detection of elements of interest by means of a 
detection system in a digital radiographic image of an object, acquired on an acquisition chain. 

[0013] According to a general characteristic of the invention, the method comprises two phases. The first is a cali- 
bration phase, in which one determines from a mathematical model of the acquisition chain and object a plurality of 
images corresponding to a plurality of gray input levels representing a set of parameters of the object and a set of 
parameters of the acquisition chain variable over a range of predetermined values. One then elaborates, from the sys- 
tem of detection and plurality of images, a mathematical model of detection giving a theoretical number of elements 
detected as a function of a set of parameters of the acquisition chain, a useful output signal of the detection system and 
a gray background level of the image considered, representing one of the sets of parameters of the object. 
[0014] Each gray level presented on input of the mathematical model of the acquisition chain and of the object 
advantageously represents a set of parameters of the object, such as the thickness and composition of a breast on a 
mammogram. 

[001 5] The mathematical model of detection thus established represents the performances of the detection system 
relative to a set of parameters variable over a range of given values. In other words, the theoretical response of the 
detection system is determined as a function of all the configurations realizable in the course of an X-ray. The first phase 
is preferably carried out only once for any detection system. 

[0016] The second phase is carried out for each acquisition of a digital radiographic image. It is a phase of use in 
which the detection system is applied on the digital radiographic image. Then, for each elementary zone of that image, 
one introduces as input data in the mathematical model of detection the useful output signal of the detection system, a 
set of parameters of the acquisition chain taken among the predetermined values and used to obtain said digital radio- 
graphic image, and a gray background level calculated on said digital radiographic image, so as to determine a theoret- 
ical number of elements detected. One then selects among the elementary zones of the digital radiographic image 
acquired those whose theoretical number of detected elements satisfies a predetermined criterion. 
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(0017J The invention will now be described in greater detail, by way of example, with reference to the drawings, in 
which: 

Figure 1 is a schematic view of the elements of the acquisition chain and of the object according to the invention; 

Figure 2 is a diagram showing the stages of the phase of calibration of the acquisition chain; and 

Figure 3 is a diagram showing the stages of the phases of use of the image. 

[0018] All of the parameters necessary for obtaining a theoretical number of elements detected for each elementary 
zone of the digital radiographic image is re-entered in the mathematical model of detection. These parameters can be 
separated into three groups, namely, the parameters of the acquisition chain, the output signal of the detector and the 
parameters calculated from the digital radiographic image, such as the gray background level. It is also possible to intro- 
duce in the parameters of the third group, in addition to the gray background level, other characteristics of the object or 
of the acquisition chain measurable on the image, like, for example, the amplitude of the periodic linear artifacts intro- 
duced in the image by an antidiffusion grid. 

[0019] The mathematical model of the acquisition chain and of the object advantageously includes a noise model. 
[0020] According to one embodiment of the invention, the detection system, adapted to each of the sets of param- 
eters of the elements of interest, is applied several times on the same image, so as to scan a given range of values of 
the parameters of said elements of interest. 

[0021 ] For example, when the elements of interest are microcalcif ications, their size and shape are preferably taken 
as parameters of those elements of interest. 

[0022] In a preferred embodiment, the detection system, adapted to each of the sets of parameters of the elements 
of interest, is applied several times on the acquired digital radiographic image, so as to obtain several theoretical num- 
bers of elements detected for each elementary zone of the digital radiographic image, and one selects among the ele- 
mentary zones the ones whose smallest theoretical number of detected elements obtained satisfies the predetermined 
criterion. 

[0023] In general, taking the smallest of the theoretical numbers is equivalent to applying a detection system best 

suited to the size and shape of the microcalcif ications. 

[0024] The detection system can advantageously be a top-hat type filter. 

[0025] The top-hat transformation is a mathematical treatment involving two stages: 

1 . Opening of the image with a structuring element; 

2. Subtraction of the image opened on the original image. 

[0026] A structuring element is a mask which, applied on an image, makes it possible to delimit a zone of pixels. 
When a structuring element is applied on an image, that is, when the mask is placed on each of the elementary zones 
of the image, all of the structures smaller than the structuring element are eliminated. A given opening corresponds to 
a structuring element of a given size. In other words, the result of the first stage is a smoothing of the gray levels. The 
gray background level is thus obtained. The second stage is an operation of subtraction of the gray background levels 
on the initial image. A high-frequency signal presenting peaks is thus obtained. This type of filtering is known to the 
expert and the latter may refer, (or more details, to the article by J. Ferra. Image Analysis and Mathematical Morphology. 
Vol. 2, Academic Press, 1988. In general, this type of filtering makes it possible to extract light and small-sized zones 
from a digital image and to clear the background. 

10027] On the high-frequency signal thus obtained, a threshold can then be applied, so as to post a theoretical 
number of elements of interest detected. Several thresholds are advantageously applied in order to cover a given range 
of threshold values. In an all-or-nothing configuration, when a threshold is applied, for every element detected, the out- 
put of the system realizing the threshold presents a signal level identical to the threshold level applied, and for no ele- 
ment detected that output presents a zero level. 

[0028] According to one embodiment, the gray background level is calculated by means of a low-pass filtering. 
[0029] According to a variant of the invention, the background level is calculated by means of a morphological 
opening. 

[0030] The theoretical number of elements detected can be advantageously expressed in terms of false positive 
probability, that is, the probability for an element to be detected when no element of interest has been introduced in the 
mathematical model of the acquisition chain and of the object. 

[0031] According to another embodiment of the invention, the theoretical number of elements detected is 
expressed in terms of false negative probability, that is, the probability for an element of interest introduced in the math- 
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ematical model of the acquisition chain and of the object not to be detected. This embodiment necessitates modifying 
the model of the object by introducing theoretical elements of interest there. 

[0032] Furthermore, satisfying the predetermined criterion involves fixing a predetermined overall threshold of 
probability and selecting elementary zones whose sum of probability values remains below or equal to that threshold. 
e In other words, the predetermined criterion can advantageously be a percentage X photon, for example. Satisfying that 
criterion boils down to preserving the elementary zones whose sum of probability values remains below or equal to X. 
The summation begins with the lowest probability values. 

[0033] Satisfying the predetermined criterion can also involve fixing a predetermined overall threshold of probability 
and selecting each elementary zone whose probability value remains below or equal to that threshold. In this case, no 

w summation is made, but all of the probability values below or equal to X are selected. 

[0034] According to an advantageous variant of the invention, the theoretical number of elements detected is 
expressed in terms of false negative and false positive probability. Satisfying the predetermined criterion involves fixing 
a predetermined overall threshold of probability and selecting each elementary zone, of which a given linear combina- 
tion of two false negative and false positive probability values remains below or equal to this threshold. 

75 [0035] Although the invention is not limited thereto, application of the method according to the invention is going to 
be described for the detection of microcalcifications by means of a top-hat transformation. 

[0036] The mathematical model of the acquisition chain and of the object is now going to be described. The acqui- 
sition chain is an image acquisition chain containing a set of elements laid out along the same axis. Referring to Figure 
1, the acquisition chain comprises an X-ray tube 1 capable of emitting a beam of X-rays 2. Placed on the path of the 
20 beam of X-rays 2 are a filter 4, an object 5 having to be studied, such as a breast, and a detector 6 connected to a digital 
image processing system 8, such as a properly equipped computer, by means of an analog-digital signal converter 7. 
This group of elements is fixed on a column (not represented) making it possible to set the distance between the tube 
1 and the detector 6. 

[0037] The tube 1 contains an anode 3 from which the X-rays 2 are emitted and a window through which the X-rays 
25 2 leave the tube. The filter 4 makes possible a bandpass frequency filtering of the X-rays 2. 

[0038] It is assumed that the breast to be X-rayed is composed of a fibrous tissue 5a (tissue a) and of a fatty tissue 
5b (tissue b). The X-rays emitted by the tube 1 are collected by the detector 6 equipped with a scintillator converting the 
X-rays 2 into visible photons and with a matrix of photodiodes in. order to obtain an electric signal. The electric signal 
thus obtained is then transformed by the converter 7 into a digital signal to be processed by the digital processing sys- 
30 tern 8. 

[0039] The X-rays 2 coming out of the tube possess an energy E. The energy is successively attenuated by the ele- 
ments crossed by the X-rays 2. These elements, the window of the tube 1 , the filter 4, the fibrous tissue 5a of the breast 
5 and the fatty tissue 5b of the breast 5 respectively possess an attenuation coefficient noted Mwindow. Miiiten tissue a and 
Ptissue b- These elements constitute a set called transmission part. When the X-rays 2 cross one of these elements, they 
35 undergo an energy loss defined by: 

* window (E) = 6*P("H window(^) ^ window) 
tfm G r( E ) s=ex P(-Mti tte r( E )- 1 fm G r) 

40 

* tissue (^,1 tissue) = ex P("M tissue 0=M tissue 

(x.y)) 

U tissue is a ,inear combination Of H tjsS ue a and ^ tissue tr 

45 t,(E) is therefore a corrective term taking into account the attenuation undergone by the X-rays 2 across the element 

i of the transmission set, 

1 1 being the distance covered by the X-rays 2 in each element i. 1 tissue is the sum of the thicknesses of the fatty and 
fibrous tissue for a point of coordinates (x,y) projected on the detector. 

50 

[0040] The energy level of the X-rays arriving on the detector 6 is therefore diminished by a quantity of energy 
according to the corrective terms t,(E) represented above and it is expressed in terms of number of X photons, noted n. 
The number n is therelore the number of X photons arriving on the detector 6 as a function of the number n 0 of X pho- 
tons emitted by the tube 1. Those numbers are functions of the energy E and of the locus of coordinates (x,y) on which 
55 the X photons arrive on the scintillator. Consequently, they are noted n(x,y,E) and n 0 (x,y.E). The number n(x,y.E) results 
from the local attenuation undergone by the n 0 (x,y,E) X-rays on their crossing between the tube 1 and the detector 6 
and depends on the terms t,(E) of each element i crossed. 

[0041] The conversion of the X photons into visible photons is made in the scintillator. The conversion of these vis- 
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ible photons into electric charges is made in the photodiodes. Then, the conversion of the charges into gray level of a 
pixel is made in the analog-digital converter 7. 

[0042] A model is used of the detector 6 having a conversion factor n{E) between the X-rays of energy E received 
and represented by n(x,y,E) and a gray level l(x,y) of a pixel of coordinates (x,y) on the detector 6. This is translated by 
the following equation: 

Emax 

I( X/ y) = J (rV5(x,y / E).n(E).twmdow(E).tfilter(E).ttissue(E / ltissue))dE 

0 

The digital resolution of that equation enables us to determine the value of thickness l t issue 
[0043] The number of X photons arriving on the detector 6 is then calculated: 

n(x,y,E) = n 0 (x I y,E).t wjnd0W (E).t filter (E).t tissue (E,1 

tissue' 

[0044] The total number N(x,y) of X-rays arriving on the detector 6 is then determined by: 

Emax 

N(x,y) =J (n(x,y,E)dE 
0 



[0045] The model of the acquisition chain and of the object therefore takes into account the acquisition chain and 
the breast, and it also takes into account the noise due to operation of the acquisition chain. An analysis of the detector 
6 makes it possible to reveal two main noise sources, the electronic noise of variance a e 2 , due to the detector 6, and 
30 the quantum noise of the X-rays 2 of variance a q 2 . 

[0046] a G is a value that is measured on the detector. The variance of several images supplied by the detector can, 
for example, be calculated when the tube 1 is not emitting any radiation. 

[0047] The distribution of the noise signal induced by the n X photons arriving on the detector 6 generally follows a 
Poisson distribution, but it is evaluated here at a Gauss distribution of variance equal to n. Therefore, a q = N(x,y). 
35 [0048] One then obtains: 

a q = vr\l(x,y). 

[0049] The noise is thus modeled by its variance o 2 : 



2 

'e 



[0050] The calibration phase is illustrated, in particular, by Figure 2. A model of the acquisition chain and of the 
object 13 is distinguished there, having as input variables a succession of gray levels 10 corresponding to different 
45 thicknesses and different compositions of the object generally expressed in terms of absorption coefficient, and the set 
of parameters of the acquisition chain, and which are variable over a range of given values: 

acquisition parameters 9: supply voltage of the tube, intensity of the supply current of the tube, exposure time, 
material of the anode, exit angle of the X-ray beam 2, 

50 

parameter of the positioner 1 1 : distance between the tube 1 and the detector 7, 

parameters of the detector 12: type of scintillator, type of photodiode. 

55 [0051 ] The model of the acquisition chain and of the object 1 3 generates several images on its output. Each output 
image can be represented by a useful signal 14. Said useful signal consists of a gray background signal corresponding 
to one of the gray input levels and of a noise signal. A top-hat transformation 15 is then applied, in which an opening 16 
makes it possible to obtain the noise signal 17 on which several thresholds 18 are successively applied. For each 
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threshold applied, the false positive probability is determined by counting the number of peaks exceeding said thresh- 
old. This operation is carried out for the different gray levels at input of the model of the acquisition chain and of the 
object so as to plot a calibration curve 19 representing the mathematical model of the top-hat transformation under dif- 
ferent conditions established by the set of parameters. This curve is in D dimensions, D being the total number of 
5 parameters considered. 

[0052] The phase of use of microcalcif ications is then undertaken in accordance with Figure 3 on a breast mammo- 
gram 20 by means of the acquisition chain 21. Several top-hat transformations 23 and 28 of different openings 24 and 
29 are applied on the image 22 obtained from the acquisition chain 21 . An opening makes possible a filtering of the sig- 
nal considered by eliminating all the elements ol size less than a given size. Openings 1 and m (m is a positive whole 
io number greater than 1) make it possible to obtain a gray background signal 25 or 30 from the image 22. A subtraction 
operation 26 or 31 is then carried out between the signal of image 22 and the gray background signal 25 or 30 in order 
to obtain the useful signal 27 or 32. The useful signal 27 is therefore the signal of image 22 minus the gray background 
signal 25. 

[0053] In each case of opening 1 to m, two signals are obtained at that stage, the useful signal 27 or 32 of image 
75 22 and the gray background signal 25 or 30 of image 22. The calibration curve 19 obtained in Figure 2 is applied. It is 
desired to obtain for each pixel, from the calibration curve actually representing theoretical values, the false positive 
(and therefore theoretical) probability, by having at the start of said curve the present acquisition parameters 33 used 
for acquisition of the digital image and those calculated from the image itself, namely, the gray background level and the 
useful signal level. The parameters to be taken into account for each pixel are a gray background level, a threshold 
20 value and the acquisition parameters: 

the value of the gray background signal 25 or 30 corresponding to the pixel considered is taken for gray background 
level; 

25 - the value of the useful signal 27 or 32 corresponding to the pixel considered is taken for threshold level; 

the parameters 33 used for acquisition of the digital radiographic image is taken for acquisition parameters. These 
parameters represent one of the sets of acquisition parameters used in the course of the calibration phase. 

30 [0054] These three groups of parameters make it possible to read on the calibration curve 19 the false positive 
probability relating to the pixel considered. The FPP (false positive probability) tables 34 and 35 relating respectively to 
openings 24 and 29 of the top-hat transformations 23 and 28 are then obtained. For at openings made, FPP tables of 
type 34 and 35 represented in Figure 3 are obtained. 

[0055] A minimization operation is then carried out in the course of stage 36, taking into account the set of FPP 
35 tables. This operation therefore consists of taking, for each pixel, the minimum value of false positive probability among 
the m values obtained. A single FPP mjn table 37 is then obtained, in which each pixel is represented by a minimum false 
positive probability value. 

[0056] The user can then intervene by indicating a number of pixels he/she wishes to keep as false positives, tak- 
ing, for example, the number 3. All of the pixels, the sum of whose false positive probability values in the FPP min table 
AO 37 is less than or equal to 3, are chosen in that case. The choice is started from the lowest values, for a low false pos- 
itive probability value signifies that there is little chance the pixel chosen corresponds to noise, that is, a high probability 
that it corresponds to a microcalcification. 

[0057] Another criterion can also be established, such as the choice of the set of pixels whose false positive prob- 
ability value is less than or equal to a given number, 0.05, for example. 
45 [0058] Furthermore, mathematical equations make it possible, in case of use of a top-hat transformation, to give 
the false positive probability as a function of a threshold s, a gray background level g and the variance of noise cr 2 . The 
false positive probability for the value of a pixel of gray background level equal to g to exceed threshold s, for a noise 
level of variance a 2 , is then obtained by executing the following algorithm: 



50 




55 [0059] 



For all of the K pixels of a given structuring element, the false positive probability is given by: 



FPP(x f y) FPP pjXGl (x,y)) K 
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[0060J The invention makes it possible to determine the microcalcrfications automatically by means of a mathemat- 
ical model of the acquisition chain and of the object, taking into account the conditions of acquisition and including a 
noise model independent of the radiographic image. 

5 Claims 

1. Method of improved detection of elements of interest by means of a detection system in a digital radiographic 
image of an object, acquired on an acquisition chain, comprising the steps of: 

io a) a calibration phase, in which: 

determining from a mathematical model of the acquisition chain and object a plurality of images corre- 
sponding to a plurality of gray input levels representing a set of parameters of the object and a set of 
parameters of the acquisition chain variable over a predetermined range of values; 

elaborating, from the system of detection and plurality of images, a mathematical model of detection giving 
a theoretical number of elements detected as a function of a set of parameters of the acquisition chain, a 
useful output signal of the detection system and a gray background level of the image considered, repre- 
senting one of the sets of parameters of the object; 

b) a phase of use in which: 

the detection system is applied on the digital radiographic image; 

25 - for each elementary zone of that image, one introduces as input data in the mathematical model of detec- 

tion the useful output signal of the detection system, a set of parameters of the acquisition chain taken 
among the predetermined values and used to obtain said digital radiographic image, and a gray back- 
ground level calculated on said digital radiographic image, so as to determine a theoretical number of ele- 
ments detected; 

30 

selecting among the elementary zones of the digital radiographic image acquired those whose theoretical 
number of detected elements satisfies a predetermined criterion. 

2. Method according to claim 1 , wherein the mathematical model of the acquisition chain and object includes a noise 
35 model. 

3. Method according to one of the foregoing claims, wherein in the course of the calibration phase, the detection sys- 
tem is applied several times on the same image, adapted to each of the sets of parameters of the elements of inter- 
est, so as to scan a given range of values of the parameters of the elements of interest. 

40 

4. Method according to any one of the foregoing claims, wherein in the course of the phase of use, the detection sys- 
tem is applied several times on the acquired digital radiographic image, adapted to each of the sets of parameters 
of the elements of interest, so as to obtain several theoretical numbers of elements detected for each elementary 
zone of said digital radiographic image, and one selects among the elementary zones the ones whose smallest the- 

45 oretical number of detected elements obtained satisfies said predetermined criterion. 

5. Method according to any one of the foregoing claims, wherein the detection system is a filter of top-hat type. 

6. Method according to any one of the foregoing claims, wherein the gray background level is calculated by means of 
so bandpass filtering. 

7. Method according to any one of the foregoing claims, wherein the gray background level is calculated by means of 
a morphological opening. 

55 8. Method according to any one of the foregoing claims, wherein the theoretical number of elements detected is 
expressed in terms of false positive probability, that is, the probability for an element to be detected when no ele- 
ment of interest has been introduced in the mathematical model of the acquisition chain and of the object. 
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). Method according to any one of claims 1 to 7, wherein the theoretical number of elements detected is expressed 
in terms of false negative probability, that is, the probability for an element of interest introduced in the mathematical 
model of the acquisition chain and of the object not to be detected. 

0. Method according to claims 8 or 9, wherein satisfying the predetermined criterion involves fixing a predetermined 
overall threshold of probability and selecting elementary zones, the sum of whose probability values remains below 
or equal to that threshold. 

1. Method according to claims 8 or 9, wherein satisfying the predetermined criterion involves fixing a predetermined 
overall threshold of probability and selecting each elementary zone, whose probability value remains below or 
equal to said threshold. 

2. Method according to one of claims 1 to 7, wherein the theoretical number of elements detected is expressed in 
terms of false negative and false positive probability, and in that satisfying the predetermined criterion involves fix- 
ing a predetermined overall threshold of probability and selecting each elementary zone, of which a given linear 
combination of two false negative and false positive probability values remains below or equal to that threshold. 



8 



EP 1 035 507 A1 



FIG. 1 
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